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lO ' Abstract 

(N ■ 

We propose an algebraic form for the density of states of 
' ch '' quarks and gluons in a Quark-Gluon Plasma (QGP)fireball 

Q_i' in quasi-equilibrium with a hadronic medium as p(k) = 

p I f + (3k + 5k 2 , and determine the parameters a, (3 and 5 

using Lattice Gauge results on the velocity of sound in 
QGP. The behaviour of the resulting p(k) can be easily 
compared with the thermodynamic data on QGP that is 
j> | expected from LHC and other RHIC experiments. Our 

ON ■ numerical result shows a linear rise of the value of p(k) for 

^ ■ k ~ T w 160 to 180 MeV, which is significant, and throws 

lO ■ light on the evolution of the QGP phase. 
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There is high expectation for early availability of data on the 
behaviour of Quark-Gluon Plasma (QGP) and its thermodynamics 
from ongoing Relativistic Heavy Ion Collisions (RHIC) experiments 
especially at LHC. Though the main theoretical results from Lattice 
Gauge Computations [1] can be tested and counter checked with 
data, there is enough opportunity to test semi-phenomenological 
models of QGP. The thermodynamics of the system can easily be 
analysed with the parametrisation proposed in the paper. Ever since 
the suggestion that the QGP can be treated as a gas in thermody- 
namic quasi-equilibrium with its hadron gas environment, there have 
been a number of proposals for modelling the gas and the density of 
states of quarks and gluons in them in the literature [3]. In natural 
units (h = c = 1), the relativistic momentum '/c' of the quarks and 
gluons is the argument of the density of states p(k), where the total 
number of states 'NqJ is given by 

POO 

n m = Kg / P( k ) dk > (!) 

J 

where the subscripts q,g refer to quarks and gluons and \ qt9 refers 
to the QCD multiplicites for quarks and gluons. we parametrise p(k) 
as 

p(k) = +f 3k + 5k 2 , (2) 
k 

where a, (3 and S are scalar parameters obviously dimensional, to 
be determined from data on QGP thermodynamics when available, 
but for the present from Lattice gauge simulation QGP behaviour. 
This structure (2), is independent of any assumption on the nature 
of the QGP droplet as in previous literature [3]. 

As in our previous modeling of the QGP droplet in the hadronic 
medium [3], we take the hadrons to be majorly made up of pions as 
the most copious component of the medium. 

The Free-energy of the system is made up of the following pieces, 
namely the free energies F q of the u,d,s quarks, the free energy F g 
of gluons and the free energy F n of the pionic hadrons. 

The total free energy F by 

F= J2 F qi + F g + F n (3) 

i=u,d,s 

for the system of QGP in quasi-equilibrium with the pionic medium. 
As in the literature [3], we have the following expression for the 
free energy pieces in (3), 
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F qi>9 = TT\ qi , 9 J dkp(k) ln(l ± e -(V^)/^) 



(4) 



and 
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= (3T/2vr 2 )z/ 



(5) 



where we take the multiplicities as usual for quark X q = 6 and 
for gluon X g = 8 and the masses m u = rrid = and m s = 150 MeV 
and v is the QGP volume which we take ~ 1 fm 3 as in previous 
paper [3] . 

With these ingredients we can compute the various thermody- 
namic parameters of the system, 

the entropy as usual is S = — (§fr)y, and the specific heat Cy = 

— The velocity of sound in the system is C 2 S = 

Thus using (4), (5) and (2) in (3) and numerically evaluating the 
integrals we have an algebraic expression of the form 



where Ct's are the numerically evaluated integrals in the expres- 
sions for S and Cy. 

In order to fix a,/3 and 5, we choose the random points in the 
C 2 vs. T simulation graph from lattice calculations [1] as shown in 
Fig. [1] and solve the three simultaneous algebraic equations using 
the values fa, fa 2 i fa 3 of C 2 S at the different temperatures Ti, T 2 
and T 3 in relation (6) of the form 



a(Ci " /t.CtO + 0(Cfi - hXl) + " MtJ = (7) 



and similar equations at temperatures T 2 and T3. We take only 
the operative region in the graph with T < 300 MeV for which 
RHIC data is expected to be available. 

Putting 




<*Cr + ffCr + $(t 



(6) 
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which can be solved for a, ft and S in terms of the known £|'s. 

To increase the accuracy, we evaluate a, ft and 5 for a number of 
random triple points in the C 2 vs. T graph with the corresponding 
fr t values in (9)and then took the average of the values so obtained. 

There is also a natural lower bound for p(k) in (2) which has 
a minimum given by a minimum energy value l k m ' satisfying the 
equation 

25k 3 m + ftk 2 m -a = (10) 

The numerical values are a = 5.7 x 10 4 , ft = —0.74 MeV~ 2 , 
5 = 1.3 x 10~ 2 MeV~ 3 and k m = 140 MeV.Fig.[2] summaries our 
result regarding the variation of p(k) vs. k dependence and the linear 
rise of the density of states for k ~ T lies in the band 160 to 180 MeV 
is significant indicating something drastic happening at these tem- 
peratures and energies. This feature can be readily checked from 
the data as soon as it is available, though this only reproduces the 
result of [1]. 
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